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Abstract 

We obtain the critical threshold for a host of Potts and percolation 
models on lattices having a structure which permits a duality consider- 
ation. The consideration generalizes the recently obtained thresholds 
of Scullard and Ziff for bond and site percolation on the martini and 
related lattices to the Potts model and to other lattices. 
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1 Introduction 



The Potts model has been in the forefront of active research for 
many years. Despite concerted efforts, however, very few exact results 
are known Unlike the Ising model for which the exact solution is 
known for all two-dimensional lattices, the relatively simple question 
of locating the critical frontier of the Potts model has been resolved 
only for the square, triangular, and honeycomb lattices |H1IU[3]. The 
determination of the Potts critical frontier for other two-dimensional 
lattices has remained very much an open problem 

In two recent papers using a star-triangle relation and a dual trans- 
formation, Scullard |7j and Ziff (Hj succeeded to determine the critical 
thresholds of site and bond percolation processes for several new two- 
dimensional lattices. As percolation problems are realized in the q = 1 
limit of the g-state Potts model [?|JE3], the new percolation results sug- 
gest the possibility that similar thresholds can also be determined for 
the Potts model. In this paper we report this extension. We derive 
more generally the exact critical frontier of the Potts model for a large 
class of two-dimensional lattices including those considered in |S] , 
and obtain the corresponding percolation thresholds. 

Consider a lattice having the structure shown in Fig. 1, where each 
shaded triangle denotes a network connected to its exterior through 3 
spins o"i,o"2,o"3. It was established by Baxter, Temperley and Ashley 
[3] using an algebraic approach that this Potts model possesses a dual- 
ity relation and a self-dual trajectory. A graphical proof of the duality 
relation was later given by Wu and Lin , and subsequently Wu and 
Zia established rigorously that in the ferromagnetic regime of the 
parameter space the critical threshold is indeed the self-dual trajec- 
tory. 

Specifically, write the Boltzmann factor for the shaded triangle as 

F(a l ,a 2 ,a 3 ) = A + B 1 S 23 + B 2 5 31 + B 3 5 12 + C 8 123 (1) 

where 6ij = S ai)a . , 5ijk = Sijdjkfiki- Then the model possesses a duality 
relation in the parameter space {A, B±, B 2 , B 3 , C} . In the ferromag- 
netic regime 

B 1 +B 2 + B 3 + C>0, B. + Bj+CyO, i^j (2) 

the critical frontier of the Potts model is given by the self-dual trajec- 
tory 

qA - C = 0. (3) 
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By realizing the shaded network as a simple triangle, for example, one 
recovers from (J5J the critical point for the Potts and bond percolation 
models on the square, triangular, and honeycomb lattices Another 
realization of the Boltzmann factor (^Q) is the random cluster model 
[Sj with 2- and 3-site interactions jllj . The isotropic version of the 
random cluster model has been analyzed very recently by Chayes and 
Lei .13. who established on a rigorous ground the duality relation 
and the self-dual trajectory ©. Our new results concern with other 
realizations of if!]). 



2 The martini lattice 



Consider the network shown in Fig. 2 as an instance of the shaded 
triangle in Fig. 1. This gives rise to the martini lattice shown in Fig. 
3 [3 |Hj . The Boltzmann factor for the network is 

g 



F(ai,a 2 ,cr 3 ) 



ex P 

{0-4,0-5, ag} = l 



Vi6u + V 2 5 25 + V 3 d; 



(4) 



+W 1 5 56 + W 2 5 46 + W 3 8 Ab + M8 A m 
where Vi nd Wi are 2-site Potts interactions and M a 3-site interaction. 

It is straightforward to cast @ in the form of (^Q) [141 115j to obtain 
A = viv 2 v 3 + v\v 2 (q + w\ + w 2 ) + v 2 v 3 (q + w 2 + w 3 ) + v 3 v\(q + w 3 + wi) 
+ (q + vi + v 2 + v 3 ) q 2 + q(wi + w 2 + w 3 ) + h 



Bi = 
C = 
where 



VjV k 



h+{q + Vi)Wi , i=£ j j 



v\v 2 v 3 h, 



(5) 



h 



e v * - 1, 

e M+W 1 +W 2 +W 3 



Wi 
,Wi _ 



li'i 



e 

*w 2 



JY 3 



+ 2. 



(6) 

As alluded to in the above, in the ferromagnetic regime Wi > 0, Vi > 
0, M > satisfying ©, the critical frontier of this Potts model is the 
self-dual trajectory (j3J) which now reads 

q(q + vi + v 2 + v 3 ) q 2 + q{w\ + w 2 + w 3 ) + /i 



+ ^1^2(^1 + w 2 + q) + ^2^3(^2 + w 3 + q) 



+v 3 vi(w 3 + w\ + q) -viv 2 v 3 h = 0. 



(7) 
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The critical frontier Q is a new result for the Potts model. 

For M = oo one retains only terms linear in h and (JJJ) reduces to 
the critical frontier q 2 + q(v\ + v 2 + v 3 ) = v\V2V 3 of the honeycomb 
lattice. For M = 0, Vi = V 2 = V 3 = V, Wi = W 2 = W 3 = W, which 
is the isotropic model with pure 2-site interactions, (0) becomes 

q{q+3v)(q 2 +3qw+3w 2 +w 3 )+qv 2 (v+6w+3q)-v 3 (3w 2 +w 3 ) = 0, (8) 

where v = e — 1, w = e w — 1. For w = v it reduces further to the 
expression 

q A + 6q 2 v + q 2 v 2 (15 + v) + qv 3 (l6 + 3u) - v 5 (3 + v) = 0. (9) 

One variation of the martini lattice is the A lattice [7J [H] shown in 
Fig. 4(a) obtained from the martini lattice by setting v% = 00, v 2 = 
i>3 = v,wi = u>2 = u>3 =to. This gives rise to the Potts critical frontier 

q 3 + q 2 (2v + 3w) + g(w 2 + 4-u;u + 3w 2 + w 3 ) - v 2 (3w 2 + w 3 ) = 0, 

A lattice, (10) 

Another variation of the martini lattice is the B lattice |Hj shown 
in Fig. 4(b) obtained from the martini lattice by setting v 2 = V3 = 
00, v 1 = v, w\ = 1V2 = W3 = w. This leads to the Potts critical frontier 

q 2 + q{v + 2w) -vw 2 (3 + w) = 0, B lattice. (11) 

Both expressions (|10|) and 1)11(1 are new. 

3 Percolation threshold 

We now specialize the above results to percolation. 

It is well-known that bond percolation is realized by taking the 
q = 1 limit of the g-state Potts model with 2-site interactions EH ■ 
For bond percolation on the martini lattice in Fig. 3, we set q = 1 and 
introduce bond occupation probabilities xi = 1 — e~ Vi , yi = 1 — e~ Wi . 
The percolation threshold then assumes the form 

^1^2(^3 + VW2 ~ 2/12/22/3) + x 2 x 3 (yi + y 2 V3 - 2/12/22/3) 
+ x 3 xi(y 2 + 2/32/1 - 2/12/22/3) - 2:13:22:3 (2/12/2 + 2/22/3 + 2/32/1 - 22/12/22/3) 
= 1 + (e M - 1)(1 - xia:2 - x 2 x 3 - x 3 x 1 + X1X2X3) . (12) 
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For isotropic bond percolation Xi = x, yi = y and M = 0, this reduces 
to the threshold 

3x 2 y(l + y-y 2 )-x 3 y 2 (3-2y) = l, (13) 

which is a result obtained in [S]. 

For bond percolation on the martini A and martini B lattices shown 
in Fig. 4, by setting y% = y and x% = x or xi = 1 (for Vi = oo) we 
obtain from l)12|) the thresholds 

2zy(l + y - y 2 ) + x 2 y(l - y) 2 = 1, A lattice, 

(l-y) 2 {l + y)-xy(2-y) = 0, B lattice. (14) 

For uniform bond occupation probability x = y = p, (|13|) and (|14|) 

reduce to 

(2p 2 - l)(p 4 - 3p 3 + 2p 2 + 1) = 0, martini lattice, 
p 5 - 4p 4 + 3p 3 + 2p 2 = 1, ^ lattice, 
(1 - 2p)(l +p-p 2 ) = 0, 5 lattice, (15) 

yielding the thresholds p c = l/y/2, 0.625457 • • • , and 1/2 respectively. 
These numbers have been reported in [71 [Hj. Note that the thresholds 
UBJ and (dU) can also be deduced from 0, (CHI), and JTTJ by setting 
q = 1, v = x /(l - x), w = y/(l - y). 

Consider next a correlated bond-site percolation process on the 
honeycomb lattice with edge occupation probabilities x±, X2, X3 and al- 
ternate site occupation probabilities s and 1. Now the site percolation 
is realized in the q = 1 limit of the g-state Potts model with multi-site 
interactions [Hlj- Therefore, by setting yi = and s = 1 — e~ M , we 
obtain from ()12|) the critical frontier for this site-bond percolation, 

S{X\X2 + X 2 X 3 + X3X1 - X1X2X3) = 1. (16) 

The expression (|16|) . which generalizes an early result due to Kondor 
|17j for xi = X2 = X3, is the central result of [7J derived from a 
star-triangle consideration. Here, it is deduced as the result of an 
application of our general formulation. 

As pointed out by Scullard [7J and Ziff the expression lfTT)|) also 
gives the threshold for site percolation on the martini lattice of Fig. 3, 
where xi,X2,x^ are occupation probabilities of the three sites around 
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a triangle and s is the occupation probability of the site at the center 
of the Y. For uniform occupation probability x\ = %2 = £3 = s, 
(|16|) yields the threshold s c = 0.764826- • • for site percolation on the 
martini lattice [7j. 

Setting X3 = 1 in (|16|) we obtain the threshold for site percolation 
on the martini A lattice of Fig. 4(a) as 

s(xi + X2) = 1, site percolation — A lattice (17) 

where x%, xi are occupation probabilities of the 3-coordinated sites and 
s the occupation probability of the 4-coordinated sites. For uniform 
occupation probability x\ = X2 = s, (|T7|) yields the threshold s c = 
1/V2 for site percolation on the A lattice. Likewise setting X2 = x% = 
1 in (|16|) . we obtain the threshold for site percolation on the martini 
B lattice of Fig. 4(b), 

s(l + x) = 1, site percolation — B lattice (18) 

where x = x\ and s are, respectively the occupation probabilities of 
the 5-coordinated sites and 3-coordinated sites. For uniform occupa- 
tion probability x = s, (fTH|) yields the threshold s c = (\/5 — l)/2 for 
site percolation on the B lattice. These results have been reported in 



4 Other lattices 



As another example of our formulation, consider the Potts model on 
the lattice in Fig. 5 with pure 2-site interactions U,V,W > 0. Writing 
u = e u — 1, v = e v — 1, w = e w — 1, we obtain after a little algebra 
the Boltzmann factor Q with 



A 
B 

C 



v 3 + 3v 2 (q + 2w) + (3d + q)(q 2 + 3qw + 3u> 2 + w 3 
uA + v 1 



2>w 2 + w 3 + (q + v)w 
u 2 (u + 3) A + 3uv 2 (u + l)(u + 2) 3w 2 + w 3 + (q + v)w 

+(u + l) 3 v 3 (3w 2 + w 3 ). (19) 

The critical frontier is again the self-dual trajectory qA — C = 0. 

The resulting self-dual trajectory assumes a simpler form for the 
percolation problem. For bond percolation we set q = 1, u = z/(l — 
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z), v = x/(l — x), w = yj (1 — y) where x, y, x are the respective bond 
occupation probabilities shown in Fig. 5. This yields the bond perco- 
lation critical threshold 



1 - 3z + z 3 - (1 - z 2 ) 3x 2 y(l + y - y 2 )(l + z) 



+ xV(3-2y)(l + 2z) 



0. 



(20) 



Setting z = in (|2()j) it reduces to the bond percolation threshold (|13|) 
of the martini lattice. Setting y = 1 l|2Ujl gives the bond percolation 
threshold 



1 - 3z + z 3 - (1 - z 2 ) 3x 2 (l + z)- x 3 (l + 2a) 







(21) 



for the dual of the martini lattice, which is the lattice in Fig. 5 with 
all small triangles shrunk into single points. 

For uniform bond percolation probabilities x = y = z = p, H2U[) 

becomes 



1 - 3p - 2p 3 + 12p 5 - 5p 6 - 15p 7 + 15p 8 



4p y = 



(22) 



yielding the threshold p c = 0.321808 • • • . Compared with the threshold 
p c = 0.707106- • • for the martini lattice, it confirms the expectation 
that percolation threshold decreases as the lattice becomes more con- 
nected. 



Summary and acknowledgment 

In summary, we have shown that the critical frontier of a host of 
Potts models with 2- and multi-site interactions on lattices having 
the structure depicted in Fig. 1 can be explicitly determined. The 
resulting critical frontier assumes the very simple form qA — C = 0, 
where A and C are parameters defined in (^Q). The corresponding 
threshold for bond and/or site percolation are next deduced by setting 
q = 1. Specializations of our formulation to the martini, the A, B, 
and other lattices are presented. 

I would like to thank R. M. Ziff for sending me a copy of [B] prior 
to publication and for a useful conversation. I am indebted to H. Y. 
Huang for assistance in the preparation of the paper. 
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Figure captions 

Fig. 1. The structure of a lattice possessing a duality relation. 
Fig. 2. The realization of Fig. 1 for the martini lattice. 
Fig. 3. The martini lattice. 

Fig. 4. (a) The martini-A lattice, (b) The martini-B lattice. 

Fig. 5. A lattice with Potts interactions U, V, W. Labels shown are 
the corresponding bond percolation probabilities x = 1 — e~ v , y = 
l-e~ w , z=l-e~ u . 
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Fig. 2 




Fig. 3 
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Fig. 5 



